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Abstract 

For the study of molecular spin junctions, we take into account two types of couplings between 
the molecule and the metal leads: (i) electron transfer that gives rise to net current in the bi- 
ased junction and (ii) energy transfer between the molecule and the leads. Using a rotating wave 
approximation in the Heisenberg representation, we derive a set of differential equations for the 
expectation values of relevant variables: electron and phonon populations and molecular polar- 
ization. A magnetic field control method to enhance the charge transfer at spin nanojunctions, 
which characterizes the molecule feature, is discussed. An approximate analytical solution of the 
resulting dynamical equation is supported by numerical solution. The magnetic control by charge 
transfer is described by transient pseudo-fermions of electrons interacting with spins. The rapid 
adiabatic passage of the energy between the molecule and the leads is taken into account. The 
current for molecular spin nanojunctions is derived. 
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I. INTRODUCTION 



A molecular spin nanojunction is a nanodevice where the electron transfer depends on the 
spin state of electrons passing through the molecule and is controlled by an external magnetic 



field. The giant magnetoresistance 



Q, 2I, y, [4] and tunneling magnetoresistance js, 



are 



widely used as molecular memory devices for magnetic field recording, and this has launched 
a new field of nanoelectronics - " spintronics" . Recently, even super magnetoresistance has 
been discussed through graphene nanoribbons|8(]. 

The main theoretical principle to describe the properties of tunneling electron transfer was 
first proposed by Gamov'Qi]. In the last few years some theoretical and computational aspects 



of electron conduction in nanojunctions have been under intense study |l 
The electric-field driven magnetic switching has been discussed 



11 



12 



HQ. 



151] and the light induced 



switching behavior in the conduction properties of molecular nanojunctions has been demon- 



strated lid. 



21 



17 



18 



19|, 



20]. The formalism of the quantized electron conductance was derived 



221 ] and the conduction behavior in the heterostructure of molecular nanojunction was 



formulated. The spin nanojunction for carbon nanotubes or graphene nanoribbons was 
investigated 23|, |2J] . Spin dependent electron transport was seen in electron nanotubes and 
graphene excited by nonstationary magnetic field |25|] as well as in graphene nanoribbons 



doped by chemically active impurities 



26 



27]. 



The formalism of the current in molecular nanojunctions was derived based on elastic 



electron scattering between two electronic baths corresponding to two leads 28 



29 



30]. 



In the elastic scattering limit, energy is lost in the lower potential lead, while no energy 



dissipation occurs in the molecular nanojunction. Following the Landauer formalism 



21 



22 



31 



321 ]. most of theoretical works on the nanojunction transport were done within 



a scattering theory approach, which neglects the contact problems and the influence of 
the scattering channels as well as the mutual influences between the electron and phonon 



subsystems [33I, l3J, l35] 



The electron transfer rate in terms of the coupling between the electronic state and the 
nuclear vibration was provided jsol. The spatial resolution at the atomic scale for single 



molecu^ 
(STM) 



es adsorbed on the surface has also been achieved by scanning tunneling microscopy 



371 ]. The nuclear vibration is described by transversal time or contact time for the 



electron transfer through a molecule. In order to estimate the transversal time, it is necessary 



2 



to take into account the chain length of a molecule in a small gap between two leads. The 



transfer time can be long enough to be comparable to the molecular vibration period 



38|. 



The strong coupling between molecular vibration and electronic states could result in trap 



ping of electrons and transition from the coherent to incoherent state 



391]. This transition 



has been found in molecular nanowires under the radiation field. It could be achieved 
by both long range electron transfers and currents through the molecular nanojunctions 
in nonequilibrium states 4^ . A class of molecules characterized by strong charge-transfer 
driven transitions into their first excited state has been investigated 28l|. The dipole mo- 
ment of such molecules changes considerably upon the excitation, resulting in a substantial 
electronic charge redistribution. 

Here, we focus our attention on the transport of the electron through the single-molecule 
spin nanojunction under the external magnetic field which controls a small nanogap between 
two metal leads. This research effort is now devoted to extending the spin-dependent ef- 
fects to magnetic molecular nanojunction for spintronic nanodevices with relatively strong 
electron-phonon coupling and large spin coherence. We investigate the spin polarized elec- 
tron transport in either occupied (OMO) or unoccupied molecular orbitals (UMO), (FigH]). 
For a molecular spin nanojunction connecting between two metal leads, time dependent 
magnetic fields can create an internal driving force for the charge to flow between the two 
leads. We suppose that the molecular junction has extremely weak spin-orbit interaction 
and weak hyper-fine interaction, meaning that the electron spin diffusion length is long 
enough to provide the spin-polarized electron injection and the spin transport between the 



leads. Our o 
and Nitzan 28 



ective in the present work is to extend the theory of Galperin, Fainberg 



29 



4l| to the case including magnetic fields and to apply the theory to the 



study of coherent control of nanojunction transport. While these problems are of general 
and fundamental interest, we note that this study is related to the efforts to develop novel 
single-electron devices, magnetic memory devices, and single-electron transistors with mag- 
netic gating[l^. In addition, the potential significance of molecular spin nanojunctions for 
device applications lies in the possibility of creating magnetic switches that could be incor- 
porated in future generations of communication systems 4^. It is conceivable that these 
devices will employ coherent spin manipulations for quantum information processing. 

In this study we particularly develop theory for the effects of strong electron-phonon 
interactions on tunnelling nanojunction and inelastic tunnelling in quantum point contacts 



3 



associated with the non-hnear conduction phenomena. There are reasons to consider molec- 
ular nanojunctions subjected to strong magnetic fields. First, the structure of such junctions 
is compatible with the configuration considered for high electromagnetic field as in tip en- 
hanced scanning near field optical microscopy (SNOM) 43|. Second, consideration of the 
spin nanojunction stability suggests that strong radiation fields should be applied in the 
sequence of well separated pulses to allow for sufficient relaxation and heat dissipation. Fi- 
nally, consideration of strong time dependent pulses makes it possible to study the way to 
optimize the desired effect for the magnetic field induced electron tunneling, i.e. to explore 
the possibility for coherent control of charge flow between the leads. 

This paper is organized as follows. SecUl] introduces our model. SecJIIII derives a closed 
set of Heisenberg equations for the average values for the magnetic spin operators and for 
the annihilation and creation operators for the electrons in the molecular states. In Sec. IIVI 
and lYl we derive the current functions for the magnetic spin operators and the formulas for 
the current and charge transfer during the magnetic pulse action, and calculate the current 
induced by a quasi-stationary magnetic pulse. The results for control of the current and the 
transferred charge by chirped pulses are summarized in Sec JVIl and the conclusion is made 
in Sec JVllI In Appendices, we show that in the absence of the radiative and nonradiative 
energy transfer couplings, the equations of motion derived here lead to the well known 
Landauer's type formula for the current. 



II. HAMILTONIAN 



As for the simplest theoretical view of an efficient spin molecular transport system, the 
actual contacts of the molecule to both electrodes are presented in Fig{Tl The metal elec- 
trodes are on the left (L) and right (R) sides. The occupied molecular orbitals (OMO) and 
unoccupied molecular orbitals (UMO) are presented by line segments, and the spins are 
depicted by arrows. A molecule is positioned between two leads represented by free elec- 
tron reservoirs L and R. The system interacts with the magnetic field. The reservoirs are 
characterized by the electronic chemical potentials fii and fiji, where the difference fii — fJ'R 
= eV is the imposed voltage bias. The coupling F is shown by double headed arrows. In 
the independent electron picture, the transition between the ground and excited molecular 
states corresponds to the transfer of an electron between OMO and UMO levels. There 
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FIG. 1: A model for the molecule-electrode junction. The electronic structure is substantially 
modified as the molecule is between two leads by a nanojunction. The right (R = \{r}) ) and left 
{L = \{l})) manifolds represent two metal leads characterized by electrochemical potentials fi^ and 
fiL, respectively. Occupied molecular orbitals (OMO) and unoccupied molecular orbitals (UMO) 
are schematically drawn. The Fermi energy of the left electrode {Efl) is higher than the Fermi 
energy of the right electrode {Efr) in the presence of magnetic field. In the Hamiltonians Eqs.(l2l 
m and[5]), indices k^k' denote the states of the metal leads, and indices m,m' denote the states of 
the molecule. Couplings and Tr are denoted by double headed arrows. 



are couplings between atomic levels and electrodes. The Fermi level of the electrodes lies 
within the gap between the highest OMO (HOMO) and the lowest UMO (LUMO) of the 
molecule. This picture assumes that the coupling between a molecule and an electrode is 
relatively weak compared with the interatomic interactions. The magnetization direction is 
either parallel or perpendicular to the external magnetic field B{r,t). In FigUl the model 
Hamiltonian Eq.([T]) consists of the unperturbed part {Hq), the perturbed part {V), and the 
magnetic field effect part (Hb)'- 

H = Ho + Hb + V, (1) 



5 



where 



'/3 



/3 



+ ^mCm_o-Cm,(T + ^kC^^^C^^a- (2) 

m&M,a ke{L,R},a 

Here, L and i? denote the left and right leads, respectively. We use h = 1 and e = 1 here 
and below. The Hamiltonian Ho (Eq.([2])) contains additive terms that correspond to the 
isolated molecule (index: m) and the free leads (index: k). Operators Cm,a (Cmo-) and Ck^a- 
(^fe cr) a^ annihilation (creation) operator of an electron in the molecule and that in 
the leads, respectively, where a is the spin (t or |); operators Oq, (a+) and 6/3 (6^) are an 
annihilation (creation) operator of a phonon in the molecule and that in the thermal bath 
or phonon reservoir, respectively. The interaction part V of Hamiltonian (Eq.([T])) is: 



V = Vm + Vm + V, 



(3) 



Here the Vm term describes the coupling between the molecular electronic subsystem and 



the free-electron reservoirs in the leads, and Vm has the form 



4^: 



vm= E (^ir^cw+/^-c.)- 

K=L,RmeM;keK,c7 



(4) 



Here h.c. denotes Hermitian conjugate. The term Vat of the Hamiltonian (EqlH]) describes 



the energy trans: 
Hamiltonian 44 , 



er between the molecule and the leads. It is written in the form of Anderson 



45 



46|: 



Vr 



N 



Yl Yl iVkk^^^4,a^k',aC+^aCm',a)- 

K=L,R m.ytm'eM;k^k'eK,cr 



(5) 



Electrons of the nanojunction are coupled to both the vibrations of the molecule and the 
electrons of the leads. The Vy term is the coupling potential of the interaction between 
electrons in the molecule and phonons in the subsystem, which is taken to be linear to the 
vibrational displacements in the form of Frohlich Hamiltonian 



46|: 



(6) 



Here, is the vibrational displacement operator: = aa + a^. Magnetic interactions in 
the Heisenberg model are employed as effective interactions between the spin of the electron 



and the magnetic fields. The model of spin Hamiltonian is described by magnetic interactions 
with a few model parameters such as spin-spin interaction {47!: 



2 ^ ^ Jm,m'SmS; 



m' 1 



(7) 



which is called Ising spin Hamiltonian [47|, |48|. Here, Sm is the localized spin magnetic 
moment, Jm,m' is the interaction potential for which each pair of spins Sm and Sm' is counted 
only once. Jm,m' is the form of Green functions consistent with the spin fluctuation theory 
based on Kohn-Sham eigenfunctionsj^o], However it is dropped out because it is a 

constant independent of the magnetic field. The Hamiltonian has the form of Zeeman 



term 



5l| with the external magnetic fie^ 



the pseudo-Fermion operators 49 



52 



53|: 



d B{r,t) for the chosen value of the g-factor with 

Hb = HB(r, t)c+^^Cm',a', (8) 

where 

HB{r,t) = g^B ^{r,t)Sm. (9) 

Here, g is the Lande g-factor which is normally close to 2; ^b is the Bohr magneton and 
i?(r, t) is the magnetic field. The c,^ „ and Cm.a operators describe pseudo-Fermion operators 
specified by a which raise the energies by ±gfiBB{r,t). Taking into account that electrons 
with the pseudo-Fermion properties exist only in the presence of the magnetic fields, we 
omit the index a in all formulas below. 



III. HEISENBERG EQUATIONS 



The physics of the system can be described by different approaches. One is the method 



of nonequilibrium Green's functions 28 



29 



54| . It has advantages in formal treatment due 



to a diagrammatic representation, and it is particularly well suited for stationary processes 
where the Dyson equation can be cast in the energy representation. For time- dependent 
processes, a method based on the Heisenberg equations of motion for the expectation values 
of the operators provides a more transparent approach, since the quantities are more directly 
related to physical observables. Such a method is adopted here. The Heisenberg equations 
for Cm(c^) and Cfc(c^) can be written as follows: 



dt 



h 



[Ho + Hb + Vm + Vn + Vv 



(10) 



I.e., 

+ ^[VW,CJ + ^[Vn,C^] + ^[Vv,Cm]. (11) 

The equation of motion derived by the mathematical transformation is presented in Appen- 
dices EHH 



K=L,Rk£K K=L,Rkjtk'&K 

and 



K=L,RkeK K=L,Rkj^k'eK 



^ / J \ ram' ~~oai"m ^ 

A. Calculation of electron transfer without external field 

We calculate the second term on the right-hand side of Eqs. f|T2!) and f|T3|) . To do this, 
we shall write down the Heisenberg equation for c^, when Vat and Vy terms are not included: 

f 4|i^o,4l4fc,41 (14) 

or 

m'GAf 

and 

k'(^K 

dc 



Since ji{c^Cm) = -^Cm + c^^, we obtain 



m'eM 
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and 



Then 



K=L,Rk'&K 



+1 V ia(*^^Vc+ c ^ 

m'&M 

4 E E vT'ifih) (19) 

where (...) denotes the averaging. Here we used {c^iCm) = {c^Cm)5mm' in the second term on 
the right-hand side of Eq. flTI?]) because of the following reasons. First, Eq. flT^ was obtained 
by neglecting Vjq when the polarization {c^,Cm)\m^m' = 0. Second, the term {c!^,Cm)\m,^m' 
is the non- resonant case. As to {c'^Ck'), this term is equal to 

(CfcCfc') = fK{£k)Skk', (20) 

where fxiE) = [exp((_E' — jj,K)/kBT) + l]"-*^ is the Fermi function, E is the energy, ks is 
the Boltzmann constant, and T is the temperature. Using (c^;C.m) = {c^Cm)5mm' = and 
Eq.(|20|) for Eq.([T9]), we obtain 

+ \vi'l''\n^-fK{eu)]. (21) 

By integrating Eq. fl2T|) . we obtain 



{ctcm) = ^Kir^ £^ dt'eM^isk - e^)it - t')][nUt') - Msk)]. (22) 
By using the formula 

f°° i ihP 

Jo " ^k 

where P denotes the principal value and by assuming that the term [rimit') — fK^^k)] is 
slowly varying as compared to the exponential function, we can move the term outside the 
integral of (Eq j2T]) . The resulting integral gives 



9 



= T^lf ''^K(t) - fK{e,)][^^ + hnSisk - ej]. (24) 

Let us show the important result of the simphfied equations for the current /. Tak- 
ing into account only Vm (Knt^^"*) Eq. llMl) . and substituting the last result into 
1 Im J2k=l,r T.k&K ^^c+c„), we have 



K=L,Rk&K 



h 

K=L,RkeK 

E [^m(^) - fKi£m)]^MK,i 
K=L,R 



where 



MK,m 



27r 



(25) 



(26) 



Here we use the formula for current /, derived from well known definitions of the molecular 
nanojunction current presented in Refs. 30|, |4l|, |55| 



^ = 4 H^'tck) = I E E(Kir^^;l.4 - h.c.) 



keL 



m=l,2 keL 

2e 



TO=i,2 fceL 



(27) 



Taking into account that {c^Ck) = {c'^Cm)*, substituting Eq. into Eq. (1271) . and using 
Eqs. (I25D and ([26]), we have 

2e 



m,=l,2 fcGL 



27r 



E yEiKir^n-"^-/-(-^M^^-^'") 



m=l,2 k€L 



J- ML,m- 



(28) 



m=l,2 
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B. Calculation of energy transfer 



We now compute the equations of motion for the relaxation induced by the molecule- 
metal lead couplings, Vm and V^. We assume that the relaxation processes due to Vm and 
Vn are not independent and also do not depend on the external magnetic field. Employing a 
Markovian approximation for the relaxation induced by the molecule-metal leads coupling, 
we derive a closed set of equations for the expectation values of binary function of the electron 
number operator in the molecule: = c+Cm or rim = {nm), where rim is the population of 
electrons in molecular state m. The polarization operator in the molecule is Bm = CmCm', 
m ^ m' and the polarization of the molecule is pm = (^a/) and = The molecular 

excitation population is hk^y = c'^/Ck- From now on, in our derivation we will consider only 
the HOMO (|1)) and LUMO (|2)) among OMO's (|m)) and UMO's (|m')), respectively, 
for simphcity, while we may introduce |1) = Y.mPMPml'm), and |2) = Y.m'GM Pm'l'm''), 



where Pm{pm') is the density of the electron level m{m') of the molecule 27|. The derivation 
could be more general so that all OMOs and UMOs might be considered. In the case 
of only HOMO and LUMO, Cm(c^) and Cm'{ c^/) will be simply replaced respectively by 
Ci(c^) and £2(0^) variables of the annihilation (creation) operators for electrons in molecular 



states 
Ref. 



1) and 1 2) [28|, |29|, |4l||. In the next formulas we shall use the notations made in 
4l|. The equations of motion in Eq. f|T2l) and Eq. fll3l) include couplings of additional 
correlations of the second order (bmk) due to the electron-transfer interaction Vm and higher- 
order correlations hMi>k,k' due to the energy transfer V^. Introducing the following notation: 
= c^ci, 6m = clc2, b^j^, = bk'k = c^Ck', hkk' = hl,^ = cl,Ck, and taking into account 
(^fcfc') = '^kk') = fk{£k)Skk', we can write: cfc^,CkC2 = cfc2C^,Ck = hM^lik- Then, Vn (Eq.©) 



is represented as 



Vat - ^ ^ {Vk!^k^,^bk'''k''bM + vl!^^}tMbk,,,k>>)- (29) 

K=L,Rk"^k"'(^K 



From the Heisenberg equation Eq.l fTOj) . we have 



dt h 

m=l,2 k"G{L,R} 



+^[VN,bMb+ ki- 



rn 
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The bilinear product of Fermion or pseudo-Fermion operators for Oj commutes with the 
bihnear product of Fermion operators for any function fj, where i j- For i 7^ j, 

[aijj] = CLifj - fjdi = 0, [nm,bMbt>k] = ^^m^A/^fc/fc - &Af&fc/fc^m = [^^m, M^fc'fc) ^ud 

[umj^jv/] = [CmCm,cfc2] = (— 1)™^"'^6m ( 1 ^ Sn^)- For the Poisson brackets for the bihn- 
ear products of the Fermion operators, we have 



K=L,Rkj^k'eK 

^ ^ X] X] ^k^kl"k"'b]i"]i"'bk',k{ ^ fim- 1). (31) 
K'=L,Rk"^k"'&K' mj^m'eM 



IV. CALCULATION OF CURRENT 



The simplest approach of transport in a molecular spin nanojunction is to assume that 
incoming electrons are scattered both at the noble metal-molecule interfaces and along the 
molecular chains. Then, the conductance will depend on the net probability of scattering 



391]. Elastic scattering does not forbid electrons to transport through the nanojunction. The 



coherent conductance takes place in most molecular chains and nanowires when the electron 
transport occurs far from a resonance frequency between the metal Fermi energy and the 
molecular eigenstates at low temperatures [391]. Landauer theory assumes that electrons 
move smoothly from one electrode (L) to another {R) only by elastic scattering within 
nanojunction. In the presence of magnetic impurities or ferromagnetic leads, electrons would 
show the spin dependent transport ( "spin valve behavior") in the presence of magnetic field. 
The spin polarized electron current emission is excited by magnetic field between metal leads. 
Thus, we take into account the linear form of new spin electron polarization added in the 
Landauer formalism. The total current I is taken by the rate of change of occupation number 



operator of electrons in the molecule is described in |30l . |41| . |55| . In Eq. (!27|) . the current / 
represents the rate of flow of electrons from the left electrode to the molecule, thus using 
the previous results for {c^Cm) we obtain 
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(it ~ 

K=L,Rkytk'eK 



(33) 

Summing Eg. (1521) and (1551) and taking the expectation values, we have: 



m'6M K=L,Rk'^K 
k'jl^keK 

K=L,Rk"^k'(^K m^m'£M keL 

for the current we have from Eg. (1271) and by full analogy with Egs. (1551) and (15^ . we obtain 
a general formula for the total current given by the rate of change of the occupation number 
operator of electrons in the molecule. 

m'€M keL 
fcGL k'y^keL 

-EE (yl^''^^ti^tci+£td2)£k')} 

keL k'^k£L 

+;^E**^(^'^)(^fc^™) 

keL 

-| E E^w(4g:c+ + c+Q^cj. 



(35) 
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Here, the energy transfer, the electron-phonon interaction as well as the external magnetic 
field effect is included. The term in the braces on the right-hand side of the last equation 
is equal to zero. This can be seen if we exchange k ^ k' in the second term of the braces. 
From Eq. fl35l) one would select the part of current excited by magnetic field by using Eqs. 
(12411281) . We have a simplified form for the current: 

Ih = —HB{r,t) 22 [^rn - fK{£k)]^ML,m- (36) 

For the Hamiltonian of magnetic field and spin-spin interaction Hs (Eq.(I71)), we have 

Is = ^ ■Jm,m'SmSm'[n'm — fxi^k)]^ ML,m- (37) 

m>m' ,acr' 

It means that there is a spontaneous current due to the transient process of spin-spin in- 
teractions with the fast damping defined by the coefficient Fml,™- The transient process of 
spontaneous spin-spin interaction current can occur in the molecular structure both in the 
presence and absence of any external field. The charge transferred during the time of an elec- 
tromagnetic pulse with the finite duration is given hj Q = I(t)dt. In the next section, 
we shall provide the computation of the electric current with {ckQaC^n) and {clQaCm)- 



V. DIFFERENTIAL EQUATIONS FOR NONEQUILIBRIUM CHARGE TRANS- 
FER 

In the form of Markovian approximation for the relaxation induced by the molecule-metal 
leads coupling for the dynamics of the electron system of the molecular junction, we have 
derived a closed set of equations for the expectation values of binary operator of {fim) = rim 
and (6m) = Pm, which are variables of the annihilation and creation operators for electrons 
in molecular states |1) and |2). Straightforward operator algebra manipulations yield and 
Pm in the rotating wave approximations. The expression for nanojunction in the Hartree- 
Fock approximation nin2 + {pm){Pm) = |(^i + ^2), is correct only when the length of 
the Bloch vector is conserved. However, due to the charge transfer between the molecular 
orbitals and metals, this value is not conserved. Thus, the so-called relevant density matrix 
of molecule pM is used as a total density matrix p, which contains the information of the 
expectation values of operators fim and pM- If chosen appropriately, the relevant density 
matrix contains the essential part of the molecular dynamics, but we employ the following 

14 



properties [53|, l56|: Tr{hih2PM) = IpmP +nin2. Due to the assumption mentioned above, 



the scattering rate of populations depends on the polarization pM which is similar to the case 



of semiconductor Bloch equations [56j]- It is assumed that the polarization is small, so that 



in the final scattering terms they keep only the terms linear to the polarization [56|. The 
advantage of our approach to the molecular junction is that it could also be used when the 
dissipative system (metals) is not in equilibrium and the many-body effects are significant. 
In this case, the total density matrix is p = pm Y\ P^^- taking into account the formal 

k&L,R 

mathematical methods, the derivation of the differential equations for the expectation values 



41 



56j. Employing Eqs.([30]), ^ and ([MD as 



from the Heisenberg equations are reported 
well as (IDip . (lD2p and (1D6P in Appendix [D] for the polarization, we finally obtain 

dpM 



dt 



-i[uo + ^ Amk - ujit)]pM 

K=L,R 

i 

I 

2 



~2^^^ ^ rim - 1) - Pm^pm 

meM 

where u(t) = uq + p(t — to), p is a. chirp rate, 

TpAf = ^ {^(rA/i^,l + TmK,2) + TNK[i^{t)]} 

1 

and 



K=L,R 



IK 



{MK)\2 



kl 



(MK)\2 



2k 



Sk - £2 

The riji is the generalized Rabi frequency |57|| consisting of two parts. 



(38) 



(39) 



(40) 



^R = gpsh ^B{t) ^ Sm 

m,m' eM 



(41) 



m>m' ,aa' 

The first part is the Larmor frequency due to nonstationary magnetic fields and the second 
part is the spin-spin interaction. The equation for the electron number has the form: 



dn. 



--nRlmpM+ ^ [/mi^(^o/2) -n^]rMif,m 

K=L,R 

+2 {rim — fK{£k)){flQm + nQm)^VK,m 

m^m'eM 

J2 {^nk{ooo)[\pm\'+ n 



K=L,R 



meM 



(42) 
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Here we introduce the next basic dynamical variables for the phonon number which is 
Uph = (a^a-Q,), and the phonon-assisted density matrices of electron-phonon interaction are 

Pa = {P^) where = i{a^ — cia). Within the Heisenberg picture, we have a set of 
differential equations for the phonon occupation number Uph and the expectation values of 



and nqm derived in refs. 



58 



59 



60 



6l|. Substituting Eqs.dCSl [Cl and [C7l ICTOl) into 



Heisenberg equations Eqs. fllOllllI) and using Pauli commutators and anticommutators Eqs. 
( lAll IA2p for the phonon-electron interaction terms, we finally obtain 

dUQm 



+2i(l — nm)Qa^QM,m 

+ {2nph + 1)(1 - 2nm){nm 



UJaUpra 

at 

2inmQa^QM,m. 



fKiek))T 



VK,m 



dt 



+2i{uJmnm - UJkfK{Bk))nQr 

+ {2nph + l)(2nm + l)(n^ - fKi£k))'^VK,r. 



dn 



Pm 



dt 



-^anqm + 2z(l - n,m)Pa^^ 



QM,m 



— fK{£k)WvK 



dnpn 
dt 



- —UJanQm + iyUJk — i^mjnpr, 
+2i{uJrnnm " ^kf K{.£k))npr, 
—2{nm — fK{£k))^VK,% 



where ujk = £k/K '^m = £m/K and 



dQa 
dt 



LOnPn 



dPa 

~dt 



-UJ, 



tQa + (^m ~ fxi^k))^ MK,'i 

K=L,R 



The equation for Up^ is 



— — fK{.^k))npm^VK,r. 



(43) 



(44) 



(45) 



(46) 

(47) 
(48) 



(49) 
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Here 



NK 



m^m'&M keK 

2n 



h 



(50) 

(51) 
(52) 

(53) 



k,k'&K 

From Eq.( |35ll we have a part of the full formula for the current taking into account the 
electron-phonon interaction: 



(54) 



VI. MAGNETIC CONTROL OF CURRENT AND TRANSFERRED CHARGE 
WITH CHIRPED PULSES 

Now we have generalized theoretical results for the charge transfer at the quasi-stationary 
strong magnetic field limit, as discussed in the previous section. The well-known procedures 
are based on the coherent excitation which produces the complete population inversion 
in an ensemble of degenerate two-level molecules with the Rabi population oscillations by 
the Gaussian pulse excitation It has been demonstrated that a molecule or an atom 



63 



64 



65 



66|. To 



excited by the Gaussian pulses behaves as a semiconductor quantum dot 
solve the main problem of the Gaussian pulse excitation of molecular levels, it requires the 
information of the resonant magnetic source, the precise control of the pulse area, and the 

ation inversion procedure, known 



chirp rate jfiTi]. In order to provide the complete popu 
as adiabatic rapid passage (ARP) [62|, |67|, l68|, l69|, [TO, ul 



72], the entire population needs 



to transform from ground |1 > to the excited |2 > electronic state. Thus, it is necessary to 
sweep the pulse frequency through a resonance. The mechanism of ARP can be explained 
by avoided crossing of dressed (adiabatic) states. In particular, starting from state |1 >, the 
system follows the adiabatic state and eventually ends up in state |2 > [t^. The scheme 
based on ARP is robust since it is insensitive to the pulse area and the precise location 
of the resonance. Therefore, we shall focus on the following ARP procedures as a way to 
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control the magnetic field induced charge transfer in molecular spin nanojunctions. Our 
formalism presented for the coherent magneto-spin properties of nanojunctions of molecular 
and quantum dots by the Gaussian pulses excitation [66], is analyzed here. 

As a particular example, we shall consider a magnetic-induced charge transfer in molec- 
ular nanojunctions. The instantaneous magnetic pulse frequency uj{t) is given as the linear 
chirped pulses uj(t) = uq — ji{t — to) (where chirped rate /i is constant) du ring the pulse 



excitation 



67l |. and the Gaussian pulse of the magnetic field is used as[6^, l67|, l73| 



B{t) = B^eM-\{5''~^l^)it-tofl (55) 

where 5 is the inverse duration of the pulse and is the amplitude of pulse. 5 and ix 
are used in dimensionless units below. Particularly interesting in this respect are molecules 
characterized by strong charge-transfer transitions that are reflected in the formation of 
an excited molecular state with a magnetic dipole. When the magnetic fields operate as 
a molecular dipole connecting two metal leads along the direction of the charge transfer 
(approximately perpendicular to the current flow axis), the magnetic pumping into the 
charge-transfer state creates an internal driving force for charge flow between the two leads. 
Here we will make a reasonable assumption that a charge-transfer transition within the 
molecule is expressed in terms of the change in relative coupling strengths of the molecular 
HOMO and LUMO to their metallic contacts. We thus investigate models in which Tmk,i < 
^MK,2, K G L, i?. This inequality reflects the fact that the excited molecular state is 
dominated by atomic orbitals of larger amplitude on one side of the molecule than on the 
other side, resulting in greater overlap with metal orbitals on that side. In the calculation 
we used the next constants 
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m- Tml,! = 0.01 eV, Tml,2 = 0.02 eV and TvK^m = 0.02 
eV. The numbers taken above for the TMK,m parameters where K = {L,R} and m = 1,2 
are reasonable, and in any case we find that similar results are obtained when they are 
changed within a reasonable range. Also, the choice 0.01 < TyK^m < 0.1 eV reflects an 
assumed lifetime of ~ 10 fs for an excited molecule at the metal surface to relax via the 
electron-phonon mechanism, which is also a reasonable number. The energy of the pulse can 



be evaluated as 1-10 eV [551]. This number is the order of magnitude of normal magnetic 
field intensities used in spectroscopy, and it should be kept in mind that it could result 
from weaker incident fields due to local field enhancement. Figs. [2lll] show the influence 



of /i or the chirp rate in the time domain on the induced current calculated by Eq.fl35 
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T 

FIG. 2: Current / (in dimensionless unit) as a function of time r = uJot for the linear chirp rate 
/i/wg = 5 • 10"3 (A), 10~2 (B), and 10"^ (C). The parameters of the calculation are as follows 
6 = 5- 10-3, huJo = 3 eV, TMR,i/nuJo = 0.04, TMR,2/nuJo = 0.03, T^/^/^o = 0.01, Tqm,^ = 0.01, 
TgA/.m = 0.99. d^Bo/hujQ = 0.2, and 5/ujq = 0.1, The picture illustrates how the current depends 
on the linear chirp rate. 

in dimensionless form: / = I{t) / eV ml,i during one magnetic pulse action. These results 
are displayed as a function of dimensionless time r = ujot with the dependence from the 
numerical solution of Eqs. fl38ti49p for a Gaussian pulse. We see that the pulse chirping 
can increase the amplitudes of the induced current, which can be explained by signatures of 
ARP (Figl2]),(Figl3l). 

The pulses of the current obtained by changing the separation of pulse compression has to 
be Lorentzian. The parameter 5 is the inverse pulse duration of the corresponding transform- 
limited pulse. The chirped frequency uj{t) changes to the resonance condition during the time 
of the pulse. Note that the local field of Eqs. ([8]) and ([9]) in the nanojunction also reflects 
mangnon excitation in the leads. The incident pulse shape affected only by the compression 
uses the possible contribution of the near-field response to both plasmonic and magnon 
excitations in the leads and to excitons in the molecule. The excitation of electric current 
by the magnetic fields in the spin nanojunction also reflects phonon-electron interaction in 
the molecule and the leads. The incident pulse shape affected only by the pulse compression 
excites the phonon-electron interaction toward the broadening of the peak of the current 
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FIG. 3: Current / (in dimensionless unit) as a function of time r = uJot for the linear chirp rate 
/i/wg = IQ-^ (A), 5 • 10^2 (B) and 10^^ (q)^ ^^ere S = IQ-"^. Other parameters are identical to 
those of Figl2j The picture illustrates how signatures of ARP increase the amplitude of the induced 
current. 



response ITSJ. Such effects are presented in Figs. [2}|11 The observable of interest is 
the magnetic field induced electronic current. Changing bias under magnetic fields with a 
fixed frequency can make the polar molecule into and out of the resonance by excitation, 
leading to highly nonlinear current voltage dependence including the possibility for negative 
differential resistance. 

Figs. E] and [H] show the calculation results of the transferred charge Q: 



as a function of the chirp rate /i in the frequency domain and 6. The calculated dependencies 
Q{fj^) of Fig l5] and Q{6) of Fig|6]are confined to the values of arguments /i and 6 corresponding 
to d^Bo/hjjjQ < 0.3 {d^ = gfiBf^ is the molecular magnetic dipole moment, Eq.Q). Since 
the theory uses the rotating wave approximation, the amplitude of magnetic field has to be 
limited. One can see that Q grows rapidly for small /i and 6. Figj5] shows that the growth 
of Q is slow for moderate 6 = 10~^ and then Q tends to become a constant value for large 





(56) 
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FIG. 4: The current / (in dimensionless unit) transferred after the completion of the pulse action 
as a function of the chirp rate fi in the frequency domain, r = LOQt. Here d^Bo/Hujo = 0.2, 6 = 10~^, 
/i/wg = 10^^ (A), 5 • 10^^ (B) and 1.5 • 10~^ (C) for the transform-limited pulse. In the course of 
chirping, the pulse energy is conserved so that J'^^B'^{t)dt = Bq is constant. Other parameters 
are identical to those of FiglSJ The picture illustrates the manifestation of ARP, i.e. how the 
amplitude of the current increases with increasing chirp rate /u. 
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FIG. 5: Charge Q (in dimensionless unit) transferred after the completion of the pulse action as a 
function of the chirp rate fi for the transform-limited pulse. Here, d^BQ/hujQ = 0.2, 5 = 10~^ (A), 
5 = 5- 10"'^ (B), and 5 = 10^^ (C). Other parameters are identical to those of Figl2j The picture 
illustrates the manifestation of ARP, i.e., how the quantity of the charge increases with increasing 
chirp rate fi. 
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0.5-/ • 
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0.00 0.02 0.04 0.06 0.08 0.10 

FIG. 6: Charge Q (in dimensionless unit) transferred after the completion of the pulse action as a 
function of 5 in the presence of energy transfer. Chirp rate is ;U = 5 • 10"'^ (A), ji = 10~^ (B), and 
fi = 5 ■ 10~^ (A). Other parameters are identical to those of FigSl 

6 = 10^^. The large pulse energy is larger than the value of Q at which the growth of 
6 slows down as presented in Figs. [5] and [6l which illustrate the influence of the energy 
transfer that diminishes the corresponding values of Q (see also Figs. [2] and [3]) • Figs. 
and [6] show the resulting behavior of the charge current induced by magnetic fields from 
Eq. fl55l) . based on the full self-consistent calculation described in Sec. II. The parameters 
used in this calculation are T = 300 K. As expected, a steady state current flows through 
the spin nanojunction in the presence of magnetic pumping. A peak of the current occurs 
at the frequency of the charge-transfer transition, i.e., the HOMO-LUMO energy gap in our 
model. The fact that magnetic current can occur in a molecular spin nanojunction with 
the postulated characteristics is a direct consequence of the fact that the charge-transfer 
properties of the molecules lead to an internal driving force that would result in magnetic 
voltage in the corresponding open circuit. Another point of concern is the thermal stability 
of spin nanojunction under the proposed thermal heating. On the other hand, the current 
calculated with these parameters (Figs|2}|l]) is of order 1 nA, implying that the magnetic 
intensity which is lower by an order of magnitude can still lead to observable currents. We 
conclude that the magnetic current in a molecular spin nanojunction is a realistic possibility. 
The values of Q shown in Fig. [6] can be rationalized by the theoretical consideration below. 
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It illustrates the influence of S of magnetic pulse, the carrier pulse frequency ui, and the 
corrected frequency of the molecular transition cjq on the transferred charge Q. As shown 
in Fig El there is an optimal parameter for magnetic pulses to provide the maximum of Q. 
This picture shows us the solution to the optimal control of the parameters of the pulse in 
order to obtain the maximum charge transfer. To end this section we note that the current 
corresponding to the expectation value of Q is 0.7 ■ 10~^^ C per pulse (corresponding to 
curve C in Figj6]) and the estimated pulse repetition frequency of 82 MHz 66|] results in a 
small but measurable value of about 10 x 10^^^ ampere. When the energy transfer between 
the molecule and electron-phonon excitations in the molecules and leads is present for the 
linear chirp, this control model can complicate the Landau- Zerner transition to a decaying 



level, which were solved in Refs. 



4l[ | and 76|. The relaxation parameters in the derived 



closed set of the equations of motion do not depend on the exciting magnetic fleld and the 
phonon numbers. The present theory could be supplied by an additional equation described 
by thermal bath, as was shown in Ref . jsi | . in the case when the Rabi frequency Qn is 
much smaller than the bath correlation frequency Uc- If molecular states Em are far from 
the Fermi levels of both leads, Uc is determined by the frequency interval for the system- 
bath interaction matrix elements V^^^^ and V^^^^ and the density of states of metal leads. 
The approximation of constant relaxation parameters, which do not depend on exciting 
magnetic radiation, is consistent with the rotating wave approximation used in our theory. 
The situation is different if we assume that the orbital energy molecular level is pinned to 
;he Fermi energy of a lead. This may lead to highly nonlinear current voltage dependence 
291]. In this case, Uc is determined also by the frequency interval at which fxi^) is essentially 
changed by ~ kBT/h. Qr can be of the same order of magnitude with Uc in the rotating 
wave approximation, and the dependence of the relaxation parameters on exciting magnetic 
fleld jsSI] has to be included in the theory. 



VII. CONCLUSION 



We have investigated a model process driven by magnetic flelds in a molecular system 
connecting metal leads. We considered the general theoretical aspects of the interaction of a 
molecular spin nanojunction with the magnetic fleld within a single-electron model. A simple 
conduction model for the HOMO and LUMO of a molecule is given in the presence of the 
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magnetic field when the electronic distribution in the molecule is far from equilibrium under 
bias voltage. This nonequilibrium state is associated with the electron flux between the 
source and leads, energy flux between the nonequilibrium molecular electronic distribution, 
and the local electronic distribution in the metals leads. This is a nonradiative dissipation 
mechanism that couples electronic excitations in the molecule to excitons in the metal. We 
have investigated two aspects of the interrelationships between electron fluxes. First, we have 
calculated the dependence of the current pulse shape of the molecule on the parameters of 
Gaussian pulse. Second, we have studied the condition under which current may be induced 
without bias voltage by the external magnetic field to adjust a nonequilibrium steady state. 
Due to the competing relaxation processes, the magnetic field induced current suggests that 
the observation is feasible. The control of electron transfer by magnetic fields in a metal- 
molecule- metal spin nanoj unction is not a simple task. The theory can be useful for the 
development of the treatment of negative differential resistance and giant magnetoresistance 



in molecular spin nanoj unctions 
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23| . The present results suggest that the experiments 



for an observation of ARP are indeed feasible. We derived the explicit solution of general 
case for the current Eq.( l35|) . It is shown that the system excited by magnetic field has an 
additional combination frequency. The differential equations Eqs. (jH]) and P6|) contain new 
additional combinational frequencies depending on the concentrations of molecular electrons 
and fermi functions. The additional combination part of the current is the effective tool to 
measure the concentration of molecular electrons and it can be extended to measuring the 
concentration of the charge carriers for another nanodevice. 

As was mentioned in SecUl future generations of the spintronic molecular systems would 
employ the coherent magnetic manipulations. We considered such coherent control pro- 
cesses for the adjustment of Q by pulse with independent parameters of fi and 6. Both 
the transfer driven current and charge that give rise to the net current in the biased spin 
nanojunction and the energy transfer between the molecule and electron-phonon excitations 
in the molecules and leads can grow with growing fi and 6 regardless of the loss of energy. It 
should be emphasized that all processes considered in this work may play important roles in 
spin nanojunction in response to the incident magnetic fields. First, direct electron-phonon 
excitations of the molecular system and the metal leads [72] may affect the response to an 
adsorbed molecule that is beyond the local field enhancement associated with the local exci- 
tation. Second, experimental realization of strong local excitations in nanojunctions requires 
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careful consideration of phonon dissipation and conduction [7a\. Phonon-electron interac- 
tion and phonon excitation may be kept under control by the junction using a sequence of 
well separated magnetic pulses, as noted in the proposed experiment. 
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APPENDIX A: HEISENBERG EQUATION FOR THE ELECTRON - 
MOLECULE INTERACTION 

We use the following Pauli commutators: 



(Al) 



[Ci, Cj] = [C+, C+] = CiCj + CjCi = 0. 



(A2) 



where and are Fermi operators, 6ij is the Kroenecker delta, and cfci = "hi, Cicf = l — Ui. 
For the unpertubed Hamiltonian Hq in Eq.([2]), we have 




(A3) 



The corresponding terms for Ck{cl^) in the Heisenberg equations is 




(A4) 



k'eK 



In Eq.f llll) the interaction term of electrons in the molecule Vm has 




(A5) 



and the interaction term of electrons in the leads has 





K=L,Rm'&M;k,k'&K 





(A6) 



m' 
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To simplify the representation of Hamiltonian Eq.(jl]) of the energy transfer, we take into 
account m G M, where M = {1,2} (where 1,2 denotes the HOMOs and LUMO). Then, the 
Poisson brackets of the electrons in the molecule have the form: 

K=L,Rk^k'£K 

= — {^kk' ''Cfc Cfc'(c^Cm + CmC^)Ci + V"^^,^ ^ ^lm{Ci Ci + CiCi)c2cliCk} 

K=L,Rk^k'(^K 

= - Yl Yl {^kk'^^^k^k'ClSim + Vi!l^ki^C2C^,ci^7) 
K=L,Rkytk'eK 

The Poisson brackets for the lead electrons have the form: 
[Vn, Ck] = [J2 Yl (^fc^^^4'Cfc'C^Ci + V^^Pctc2C^,Ck"), c^] 

K=L,Rk"jtk'f^K 

= Y Y {^fc'^^Hcfc"4'C^CiCfc - CfcC+,,Cfc/c^Ci) + V;^,^5^(c 

K=L,Rk"j^k'€K 
- Ckcfc2cl,Ck")} 

= Y, Y^ {V^!^^f\cpCk'Ck - CkC^,,Ck')c2Ci + V^!lPciC 
K=L,Rk"^k'&K 

-{ Y ^kk'^^^k'ctci+ Y ^kk'^^^tc2Ck"} 
k'ytkeK k"^k&K 

= - Y ^k^^\^k'ctci+ctc2Ck') = - Y ^kk^^\^tci+ctc2)Ck'. 
k'^kf^K k'^k&K 

To obtain the formulas presented above, we used the following rules of the triadic multipli- 
cation for operators. If A; 7^ A;', then {c'^„Ck'Ck — CkC'l„Ck') = —{cl„Ck + CkC'l„)ck' = —5kk"Ck'- 
\i k = k', then {clnCkCk — CkC'^„Ck) = —CkclnCk = cl„CkCk = 0, since k" ^ k'. Therefore, we 
have {cpCk'Ck - CkC^Ck') = -5kk"h'k"Ck' for k" ^ k' . If k ^ k", then {c^,Ck"Ck - CkC^,Ck") = 
-{ctfCk + Ckcl,)ck" = -Skk'Ck"- If k = k", then {cyCkCk- CkCyCk) = -CkC^,Ck = 0. Therefore, 
we have {c^,Ck"Ck - CkC^,Ck") = -6kk'Skk"Ck" ■ 
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APPENDIX B: HEISENBERG EQUATION FOR THE ELECTRON - MAG- 
NETIC FIELD INTERACTION 



In the Heisenberg equation with the magnetic fields (llip . the term of Hb in Eq.([H]) for 
the interaction of the electron in the molecular system has the form: 

Hb = -HB{r,t){c^Ci + c^Cs) 
= -^{9f^B ■ e{c2CiB{t) exp[-iujt + i(f{t)] 

+ctc2B*{t) exp[iujt - i^it)]} (Bl) 

where e is the unit vector. Substituting Eq. flBip in the Poisson brackets of Eq. ffTT]) . we have 



[Hb, Cm] = -^{gf^B ■ e)[ctciB{t) exp[-iujt + iip{t)] + ctc2B*{t) exp[iujt - iip{t)], Cm] 

= -^{9f^B ■ e){(c^cic„ - CmC2Ci)B(t) exp[-iujt + iipit)] + 



[clc2Cm - CmCiC2)B*{t) ex^[iujt - I'^it)]} . (B2) 



Taking into account that m G M, M = {1, 2} for Eq. (lB2[) . we have: 

[Hb, ci] = ^(fi'/ifi ■ (i)B*{t) exp[iu;t - i'^{t)]c2 

and 

[Hb, C2] = ]^{g^i-B ■ e)B(t) exp[-iLjt + i(pit)]ci. 

Similar to this, the Poisson brackets of Eq. (ITTl) for the magnetic field and the lead electrons 
yield: 



[HB,Ck] = -Hij(r,t)[c^ci + c+C2,Cfc] 
-HB(r,t){(c+ciCfc - Ckcjci) + (c^CaCfe - CfeC^Cs)} 

-HB(r,t){(c^Ci - C^Ci)Ck + (C|C2 - C^C2)Ck} = 0. 
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APPENDIX C: HEISENBERG EQUATION FOR THE ELECTRON - PHONON 
INTERACTIONS IN THE MOLECULE 



In the Heisenberg equation fllip . the term of Vy in the interaction of molecular phonons 
and electrons Ck has the form: 



dcr, 



dt h 



dt h 



Vy, Cm 



(CI) 



(C2) 



Following the phonon commutation rules [ao,^^] = 1 for the Poisson brackets of the inter- 
action of molecular phonons and electrons, we have 



(C3) 



meM,k<=K 

meM,keK 



(C4) 
(C5) 



meM,keK 

+ E 2{emnm- £kfKi£k))c'lP^Cm 
■m£M,keK 



(C6) 



K=L,RmeM,k€K 



-(1 - fK{£k))nm)Qa, 



(C7) 



2 E E V^r\{l-nm)fK{ek) 

K=L,Rm£M,kGK 

-(1 - fK{ek))nm)Pa, 



(C8) 
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VM,ctP:Cm 







2 



2i V V^^'^^c+c 



(C9) 
(CIO) 

(Cll) 



(C12) 



APPENDIX D: CALCULATION OF (cjcm) IN THE PRESENCE OF A MAG 



NETIC FIELD 



In this section we calculate the expression for the expectation value (c^Cm) by taking 
into account the magnetic field, namely using Eqs. (IT^ . (IT^ . (155]) and For m = 1, we 
obtain 



m'=l,2 K=L,Rk'£K 

+ ^ E V^r\4c^a2Ci) + ^(^?/.B-e)S*(t)expM-.V^(t)](c+C2) (Dl) 

~i E E (^ii^^C+C2C+Cfc"). 
K=L,Rk"^k'eK 

Since the term (c^/Cm) 1^^^^' is the non-resonant case for multiplication of ek — S2 and e^ — £i, 
we must put {c'^,Cm) = {c^Cm')Smm' in the second term on the right-hand side of Eqs. flDip 
and dm]). For m = 2, 
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m'=l,2 

4 E E i^S^'i^^s^') 

(D2) 

In addition, using Eq.(l20l) and disregarding terms of ~ \/(^^) as a first step, we obtain 

+ '^(9I^B ■ e)B*{t)exp[iujt - iip{t)]{c'lcm'^m) (D3) 

In order to derive the solution of Eqs. flDlllD2p we shall use the slowly varying amplitude 
method. Lets us define 

{ctcm) = {c'lcm)oexp[iujt - i^it)], (D4) 
by substituting Eq. flD4p for Eq. flDSp . we have 



^(^t^ra) = ji^k- ei)(c+c^)oexp[ia;t - i^{t)] + ^V^l^^n^ - ^^^^^^ fK{ek) + 

Conserving only the resonance terms in Eq.( ]D5[) . we have 
d i 

^(c^c„)o = ilisk - e„,)/fi- uj{t)\{clcm)Q + ^idf^B ■ e)B*{t){clcm'^ra)- (D6) 
Integrating the last equation, we obtain 



{clc^)o = —{g^iB- fi)B*{t){clcra'^ra) j rfr exp ( (^fc - - )r] , (D7) 
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and 



(c^c^) ^ jV^l''\n^{t) - fK{ek)\ I drexp[-{e, - 5^)r]- (D8) 

- ■ e)^|S(t)p(c+c„) / cirexp[-(£fc -e„)r] / dTexp[i{{ek - e^'jc^)/h- uj{t))T] 

(D9) 

or 



1 i f°° 

{I - ^{gfiB ■ ef\B{t)\^ dTexp[-{ek- em)r] / dr exp[i{{ek - em'^m)/fi - i^{t))r]} 

(DIO) 

The last equation gives "saturation" corrections to the prehminary result Eg. (1241) by tak- 
ing into account the external magnetic fields. The approximation for the magnetic field 
dependent series could be produced by formal integration of Eq. ( IDIOI) with an appropriate 
treatment related to the magnetic field dependent relaxation parameters, namely 

{ctCm) = Mt) - fK{ek)){l - mt)\^TBK,mWMK,m, (Dll) 

where 

J- BK,m — 7 J- MK,m- \^^^) 

APPENDIX E: SIMPLIFIED EQUATIONS OF CURRENT FOR THE Vm IN- 
TERACTION 

In this section we derive the simplified formula for the current and set the differential 
equations for rii and n2 without Vy, Vn and He- Using the formula for the current Eq. (!35!l . 
we obtain 

/ = e{[-fL{ei) + n,]TML,i + [n2 - fL{e2)]TML,2}- (El) 
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For K & L,R, the damping coefficient has the form: 



Tmk,! = ^ E - ^o/2), (E2) 



h 

keK 



Tmk,2 = ^ E \^2T^'^^^^k - ^o/2). (E3) 



h 

keK 



The ordinary differential equation Eq. (H2l) has the simphfied forms: 

drii 



dt 



VtR\mpM+ ^ [fiKi^o/'^) -ni]TMK,i (E4) 



K=L,R 

^ = nRlmpM- E [^2 - /2x(^o/2)]rMK,2- (E5) 

K=L,R 

The differential equation for the polarization Eq. fl38|) has the form: 

^^ = -z[cJo+ AMK-^^it)]pM-^^Rini+n2-l)-^ {'rMK,2 + '^MK,l)PM- (E6) 

K=L,R K=L,R 
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